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Abstract. Barchan dunes are found where sand availability is low and wind direction quite constant. The 
two dimensional shear stress of the wind field and the sand movement by saltation and avalanches over a 
barchan dune are simulated. The resulting final shape is compared to the results of a model with a one 
dimensional shear stress. A characteristic edge at the center of the windward side is discovered which is 
also observed for big barchans. Diffusion effects reduce this effect for small dunes. 

PACS. 47.54. +r Pattern selection and formation - 45.70. -n Granular systems - 92.10. Wa Sediment trans- 
port - 92.60.Gn Winds and their effects - 91.10.Jf Topography; geometric observations - 92.40.Gc Erosion 
and sedimentation 



1 Introduction 

Barchan dunes are highly mobile dunes which constitute a 
considerable threat to infrastructure in arid regions with 
sand. They move over roads, pipelines and cover even 
cities. The large time scales involved considerably com- 
plicate reliable measurements. Nevertheless, over many 
decades measurements have been made of dunes all over 

In order to learn more about the mechanism driving dune 
morphology and dynamics numerical models have been 

els revealed many interesting results and rised new ques- 
tions which have to be answered. The inclusion of the lat- 
eral component of the shear stress of the wind field over a 
barchan dune leads to the results presented in this paper. 

The following section will explain some aspects con- 
cerning barchan dunes. The model will be introduced in 
order to simulate barchan dunes. In the following section 
scaling relations calculated with the model will be pre- 
sented and compared to measurements. The difference be- 
tween the numerical results and real dune shapes for small 
barchans will lead to an inclusion of diffusion in the salta- 
tion transport calculation. The last section will be about 
barchanoids. 



2 Basics 

The word barchan comes from the turkish language and 
means "active dune" . It was preserved in the scientific lit- 
erature to name the isolated crescent-shaped mobile dune. 
Less than 1% of all dune sand on Earth is contained in 
barchan dunes. These dunes exist mainly in areas where 
not very much sand is available and wind stays unidirec- 
tional. 



The size of barchans varies from heights of some me- 
ters (Figure ^1 to over 50 m (Figure EJ . Barchan dunes 
are not exactly shape invariant. There is a minimal height 
of 1-2 meters below which it looses its sand and no sta- 
ble shape is reached. Dune shapes seem to be controlled 
by the saturation length of the saltation transport on the 
windward side |23liU] . This is the distance needed to reach 
the saturated sand flux over a sandy surface. The satura- 
tion length, not to be confused with the saltation length 
(mean length of a grain trajectory in air), has a complex 
dependency on the air shear stress. Zero flux over a ground 
without sand bed needs a transient length to reeach a sat- 
urated flux at the place where the surface begins to be 
covered with sand. 

Small stable barchan dunes have a short slip face and 
the crest does not coincide with the brink. Whereas the 
sand is trapped completely in the slip face, the lack of a 
slip face at the horns allows the sand grains to leave the 
dune there. Thus a barchan dune grows if sand influx is 
larger than outflux and vice versa. 



3 The model 

The model described here can be seen as a minimal model 
including the main processes of dune morphology. As pre- 
decessors of the model described here the works of 21 and 
|25j revealed interesting new insights into dynamics and 
formation of dunes. The model of |2S], modelling trans- 
verse dunes, will be used here to obtain the steady state 
of barchan dunes. The simulations are carried out with a 
completely unidirectional and constant wind source. In 
this article the model is used to simulate isolated sin- 
gle barchan dunes. In every iteration the horizontal shear 
stress T of the wind, the saltation flux q and the flux due 
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ary layer as a basis of our shear stress calculation, 



Fig. 1. Aerial photography of the dune field north of 
Laayoune, Morocco. The dune heights are of some meters. The 
wind is blowing from NNE to SSW (photo taken from Sauer- 
mann (2001)). 




Fig. 2. Aerial photography of the dune field of Jericoacoara 
near Fortaleza, Brazil. Some dunes are over 50 meters high. 
Photo taken from Jimenez et al (1999). 



to avalanches are calculated. The time scale of these pro- 
cesses is much shorter than the time scale of changes in the 
dune surface so that they are treated to be instantaneous. 
In the following the different steps at every iteration are 
explained. 



v(z) ~ — In — , 

K Zo 



(1) 



where v{z) is the horizontal wind velocity, zq the rough- 
ness length, giving a measure of the roughness of the sur- 
face, and — \Jt I p is called shear velocity. Although 
it has the dimension of a velocity the shear velocity m* is 
used as a measure for the shear stress. The shear stress 
perturbation over a single dune or over a dune field given 
by the height h{x^ y) is calculated using the algorithm of 
|32j . The Ta;-component points in wind direction and the 
Tj,-component denotes the lateral direction perpendicular 
to the wind. The calculation is made in Fourier space, kx 
and ky are the wave numbers, 
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where |fc| = ^k^ + k^ and 7 = 0.577216 (Euler's con- 
stant). U{1) is the normahzed velocity of the undisturbed 
logarithmic profile at the height of the inner region I |5J 
defined in The values of the roughness length zo and 
the so called characteristic length L were adjusted to the 
data of |211 supposing we have an averaged shear velocity 
of u, = 0.5 ms~^ in the region where the measurements 
were carried out. Thus zq is set to 0.0025 m and L to 10 
m. 

Equations Q and © are calculated in Fourier space 
and have to be multiplied with the logarithmic velocity 
profile of Equation in real space in order to obtain 
the total shear stress. The surface is assumed to be in- 
stantaneously rigid and the effect of sediment transport is 
incorporated in the roughness length zq. For slices in wind 
direction the separation streamlines in the lee zone of the 
dunes are fitted by a polynomial of third order attaching 
to the brink continuously. The length of the separation 
streamlines is determined by allowing a maximum slope 
of 14° 125. The separation bubble guarantees a smooth 
surface for the wind field calculation and the shear stress 
in the area inside the separation bubble is set equal to 
zero. Problems can arise due to numerical fluctuations in 
the value of the windward slope at the brink where the 
separation bubble begins and its influence on the calcu- 
lation of the separation streamline for each slice. To get 
rid of this numerical error the surface is Fourier-filtered 
by cutting the small frequencies. 



The air shear stress t at the ground: We use the well The saltation flux g: The product of the sand density and 
known logarithmic velocity profile of the atmospheric bound- the sand velocity leads to the sand flux over the surface, 
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q(a;,y) = u{x , y) p{x , y) . The saltation flux is calculated where Vcs is the velocity of the grains in the saturated 
from mass conservation, state, 

dp(x,y,t) „ , , , , 2u, 



+Cd^ff^p{x, y, t) = r{x, y, t), 



(4) 



where p is the density and u the velocity of the sand grains 
in the saltation layer. These variables are integrated over 
the vertical component. Cdi// denotes a diffusion constant 
which will be the crucial parameter in Section [S] In the 
other chapters Cdiff is set to zero. Grains are entrained by 
the air and deposited on the sand bed which is expressed 
by the exchange term 7^. The velocity of the sand grains 
is calculated from momentum conservation, 



d\i{x,y,t) 
dt 



{\i{x,y,t)V)u{x,y,t) 




and 



(9) 



(10) 



In Equation (jS]) the force terms are arranged in the same 
order as in Equation (jSJ . The constants and model param- 
eters have been taken from 23 and are summarized here: 

g = 9.81mS"2,K = 0.4, Pair = 1.225 kg m"^, Pquartz = 

2650kgm-3, = 0.04 m, zq = 2.5 10^^ m, D = d = 
250 Aim, Cd = 3, u^t = 0.28ms-\ 7 = 0.4, a = 0.35 and 
zi = 0.005 m. 



p{x,y,t) 



{{dragix, y, t) + ibed{x, y, t) + ig{x, y)) . (5) 



The drag force ^drag denotes the interaction between air 
flow and sand grains, fbed is the decceleration of the grains 
by the collition with the sand bed and ig is the grav- 
ity force which drives the grains into the direction of the 
steepest gradient of the height profile. The time to reach 
the steady state of sand flux over a new surface is sev- 
eral orders of magnitude shorter than the time scale of 
the surface evolution. Hence, the steady state is assumed 
to be reached instantaneously. The characteristic length 
scale to reach a saturated saltation layer influences the 
calculation by breaking the scale invariance of dunes and 
by determining the minimal size of a barchan dune ^23| . 

A calculation of the saltation transport by the well 
known flux relations [Tlll7ll27j would restrict the model 
to saturated sand flux which is not the case observed at 
the foot of the windward side of a barchan dune due to 
the reduced sand supply. We simplified the closed model of 
pn] , resulting from Equation and ©, by neglecting the 
time dependent terms and the convective term of the grain 
velocity u(a;,y). This yields Equations © and ^ where 
p and u are determined from the shear stress obtained 
before and the gradient of the actual surface, 



div(pu)-,Q.,,^,.lp(l-^) jf'^) ^<^^ 

Ts \ Ps / [ 1 P> Ps 



with 
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Ps = — (rl-T-t) Ts = J— -. (7) 
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and 

4 Pquartz 
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Avalanches: Surfaces with slopes which exceed the max- 
imal stable angle of a sand surface, the called angle of 
repose O ~ 34°, produce avalanches which slide down in 
the direction of the steepest descent. The unstable sur- 
face relaxes to a somewhat smaller angle. For the study 
of dune formation two global properties are of interest. 
These are the sand transport downhill due to gravity and 
the maintenance of the angle of repose. To determine the 
new surface after the relaxation by avalanches the model 
proposed by |2| is used. For the equations see Like 
in the calculation of the sand flux the steady state of the 
avalanche model is assumed to be reached instantaneously. 
Hence, we can neglect the time dependent terms. In the 
dune model a certain amount of sand is transported over 
the brink to the slip face and in every iteration the sand 
excess is relaxed over the slip face by this avalanche model 
determining the steady state. 



The time evolution of the surface The calculation of 
the sand flux over a not stationary dune surface leads 
to changes by erosion and deposition of sand grains. The 
change of the surface profile can be calculated using the 
conservation of mass. 



dh _ 1 

dt Psand 



Vq. 



(11) 



(8) 



Finally, it is noted that Equation is the only remain- 
ing time dependent equation and thus defines the charac- 
teristic time scale of the model which is normally between 
3-5 hours for every iteration. 



The initial surface and boundary conditions: The gen- 
esis of a dune is still not known very well. Hence, the 
simulations must start with an essentially arbitrary ini- 
tial surface. We can then observe how the system reaches 
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a final steady state. A steady state means that the dune 
shape does not undergo temporal changes anymore. 

As initial surface we take a flat bed of solid ground 
with a Gaussian hill of sand on it. The maximum slope 
of the Gaussian is restricted in order to stay within the 
approximations of the shear stress calculation. 

We use quasi-periodic boundary conditions, i.e. the 
sand outflux at the outlet of the simulation field is inte- 
grated and the same amount of sand enters the field inlet 
as a constant influx. Neglecting the amount of sand leav- 
ing the simulation area by the lateral boundary, the sand 
volume is conserved. 



4 Scaling laws 

In this section the morphologic relationships between height 
h, width w, length / and velocity Vd of barchan dunes re- 
sulting of the simulations are presented. The shapes for 
different dune sizes are compared. Therefore calculations 
of dunes which have different sizes were performed for 
the shear velocities = 0.4 ms~^, u^, = 0.45 ms~^ and 
It* = 0.5 ms~^. Single isolated barchan dunes are mod- 
eled using a quasi-periodic boundary until the final steady 
shape is reached. 



Height, width and length relationships: Linear relation- 
ships between height, width and length were observed by 
0] and for barchans in southern Peru and by [S] in 
Morocco. Figures 121 and 01 depict the height width and the 
height length relationships for different shear velocities, 
respectively. 

A linear relationship is obtained for dunes larger than 
2-3 meters. The linearity is valid only for length scales 
much larger than the saturation length (I 3> Is) - The func- 
tional dependence of the saturation length on the shear 
velocity predicts larger width and length for a dune of 
the same height, for decreasing shear velocities. This is in 
agreement with !J3; who found an increase of Is at the 
shear velocity threshold. 

The asymptotic slopes of the linear relationships are 
not affected by different shear velocities. The results are 
compared to the data of [5| {Figures|21and^. The data of 
the barchans of Morocco fit best to the simulation results 
for a shear velocity u^, = 0.5 ms~^. Table ^ compares the 
slope aw and the axis intercept bw of the linear relation- 
ships to the data of [H|, [HI and @]. 
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Fig. 3. Height width relation of dunes in Morocco (diamonds) 
and the results of numerical calculations for different shear 
velocities. 



where $dune is in principle dependent on the shear ve- 
locity. The relation between dune velocity and height re- 
vealed the same deviations from Bagnold's law for small 
dunes as in the case of the height length relationship (Fig- 
ure ^J. Figure [3 depicts the dune velocity Vd versus dune 
length /. The surprising result is that the dune velocity 
does not differ from Equation Eleven for small dune sizes 
and so the saturation length seems to have no influence 
on the dune velocity. 



The shape: Barchans are not really shape invariant due 
to the important role of the saturation length of the salta- 
tion transport. Nevertheless the shapes of barchans of dif- 
ferent sizes are compared in order to obtain more informa- 
tion about the deviation from shape invariance. !5: fitted 
the normalized transversal profile of the measured dunes 
in Morocco with a parabola. 21 found a good agree- 
ment of the normalized transversal profile of his numeri- 
cal calculations (neglecting the lateral shear stress) to a 
parabola. Following his work the axes are rescaled and 
dimensionless variables are introduced, 



1 



1 



1 



x = -X y=—y z= -z. 
L w n 



(13) 



Dune velocity: According to [23 the dune velocity de- 
creases inversely proportional to the length of the enve- 
lope of the surface formed by the height profile and the 
separation bubble. The separation bubble of the barchans 
in the dune model described here fills almost exactly the 
region between the horns so that the length I of the dune 
can be used in order to evaluate the simulation results. 
Thus Bagnold's law, the reciprocal proportionality of the 
dune velocity Vd to the height h, has to be modified, to 



Table 1. Height and width relationship. 



(12) 





aw 


bw 


Finkel 


10.3 


4.0 


Hastenrath 


8.2 


9.5 


Herrmann 


11.1 


5.6 


Simulations: 






u* = 0.4 ms~^ 


9.9 


16.6 


u, = 0.45 ms"^ 


10.6 


10.4 


Uf = 0.5 ms"'^ 


10.5 


8.5 
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Fig. 4. Height length relation of dunes in Morocco (diamonds) 
and the results of numerical calculations for different shear 
velocities. 




Fig. 5. The velocity Vd of barchans fits very well to their recip- 
rocal length. Note: The velocity is assumed to be constant for 
365 days of wind per year and is thus smaller in real conditions. 



Figure El shows our normalized transversal profiles com- 
pared with the field data of PI]- In Figure [7| we compare 
the results of the normalized transversal profiles with a fit 
to a parabola. The upper part of the profile fits quite well 
whereas the lower part is far away from a parabola. We 
searched for a better fit with another function. Figures 
and IHl indicate that the powers of a cosh-function fit quite 
well, even the longitudinal profile. The inclusion of lateral 
shear stress in the model seems to lead qualitatively at 
least for the transversal cuts to a change from a parabolic 
to a cosh^-profile. In Figure |H1 the slope at the brink at 
the windward side decreases for larger dunes. Finally the 
shapes of a 45 and a 5 meter high dune are depicted in 
Figure The slip face of larger dunes cuts a larger piece 
from the dune body. 

In figure El we can see a rather sharp edge in the shape 
at the center of the windward side which is not observed 
in the results of the dune model of [23 ■ The inclusion of 




0.5 1 1.5 2 2.5 

normalized width y 
Fig. 6. Normalized superposed transversal profiles of the sim- 
ulations (solid lines) and normalized transversal profiles of the 
dunes measured in Morocco from Herrmann and Sauermann 
(2000) (dashed lines). The shear velocity in the simulations is 
It, = 0.5ms~^ . 




0.5 1 1.5 

normalized width y 
Fig. 7. Normalized transversal profiles of the dunes. A fit with 
cosh^ (x) is given by the dashed line and fits much better than 
a parabola (dotted line). The shear velocity is u, — 0.5ms^^. 



the lateral component in the shear stress calculation seems 
responsible for this difference. This edge has not been ob- 
served for small dunes, for example in Morocco. But the 
large dunes in Figure (21 show a very similar edge. Exact 
measurements of the shape of a large barchan would reveal 
more information and hopefully validate the model results. 
We assume that the absence of the edge for small dunes 
can be explained by diffusion-like smoothening effects oc- 
curing during the saltation transport. Diffusion which in 
the simulations depicted so far has not been considered 
acts on small scales and for larger scales the influence of 
diffusion on the dune shape should be negligible. In the 
next section the shapes obtained from the model are stud- 
ied for different diffusion constants. 
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0.5 1 1.5 

normalized length x 
Fig. 8. Normalized longitudinal profiles of the dunes. A fit with 
cosh*(a;) (dashed line) reproduces quite well the windward side. 
The shear velocity is Ut = 0.5ms~^. 



5 The effect of diffusion 

The dune model including the lateral component of the 
shear stress revealed a characteristic shape with a rather 
sharp edge on the windward side which is not found for 
small barchan dunes for example in Morocco. To better 
understand this effect we simulate barchan dunes of about 
12 meters height having the same volume with different 
diffusion constants Cdiff- For larger dunes and therefore 
larger length scales the effect of diffusion transport will be 
weaker and the final state calculated by the dune model 
without diffusion should match real dune shapes. 

The results in this section can be discussed only qual- 
itatively due to the fact that experimentally the diffusion 
constant is unknown. Figures 1101 and 1111 show the con- 
tour lines of the height profile and the sand density in 
the saltation layer, respectively. Very high diffusion con- 
stants like the one of the bottom part of these figures show 
that the sharp edge in the shape of the windward side is 
smoothened. The sand density in the saltation layer is flat- 
tened out as can be expected from diffusion processes. Be- 
hind the brink always the sand density drops to zero due 
to the vanishing shear stress. Diffusion changes the con- 
tour lines of the shape dramatically. An estimation of the 
diffusion constant and precise measurements of the height 
profiles of small and large barchans could give more in- 
sight into the infiuence of diffusion on the morphology of 
barchan dunes in the future. 

Figures IT?1 and lTSI depict a longitudinal and a transver- 
sal profile of a barchan dune in steady state for different 
diffusion constants. The height and the slope at the brink 
decrease for a stronger diffusion. The length of the lon- 
gitudinal profile increases for a larger diffusion constant 
whereas the width stays constant. The knowledge of Cdiff 
would also help to determine more realistic values for zq 
and L (Equation (0)) by fitting them to simulation re- 
sults from a model with diffusion. Consistency can be 
checked by comparing with higher barchans where diffu- 
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Fig. 9. On top the surface of a 45 m high barchan dune is 
depicted. The dune in the bottom part has a height of 5 m. 
There is no simple shape invariance. The shear velocity is it, — 
0.5ms~^ for both simulations. 



sion is weaker. But the computional costs and the missing 
measurements for large barchans (there do not exist many 
large barchans) make the realization still difficult at this 
point. 
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Fig. 10. Comparing the topography of two barchan dunes of 
equal volume in steady state calculated from the same ini- 
tial surface for a diffusion constant of Cdiff ~ 0.0 (top) and 
Cdiff = 4.0 m^s~^ (bottom). The shear velocity is u, — 
0.5ms-\ 
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Fig. 11. Comparing the sand density p of the saltation layer of 
two barchan dunes of equal volume in steady state calculated 
from the same initial surface for a diffusion constant of Cdi ff~ 
0.0 (top) and Cdiff = 4.0 m^s~^ (bottom). The diffusion term 
smoothes the sand density in the saltation layer. The shear 
velocity is u* — 0.5ms~^. 



6 Barchanoids, between barchan and 
transverse dunes 

If an area does not contain sufficient sand to form trans- 
verse dunes but too much sand is available to keep the 
barchans isolated from each other the barchan dunes con- 
nect and interesting hybrid forms appear, the so called 
barchanoids. The barchan dunes can be connected longi- 
tudinally and laterally. 

Here only qualitative results will be presented. The 
simulation with the dune model is performed with a nearly 
flat initial surface (some small Gaussian hills randomly 
placed on a plane) of an averaged sand height of three me- 
ters. The boundary conditions are quasi periodic in wind 
direction and open in the lateral direction. The volume of 
sand in the simulation is held constant. Figures ITU and [T51 
show two states of the simulation. In the first figure the 
barchans are connected longitudinally and laterally. The 
barchans are growing so that the sand finally accumulates 
in a single row of larger dunes which are connected later- 



ally. Similar barchanoids can be found for example in the 
dune field of Lengois Maranhenses in Brazil where on top 
in the rain season the dunes can be separated by lagoons 
filled with rainwater (Figure EJ). In this dune field the 
barchan dunes are connected in both directions. 



7 Conclusions 

The scaling laws revealed approximately linear relation- 
ships between height, width and length of the dunes. This 
is not valid for small dunes due to the influence of the sat- 
uration length. Bagnold's law was modified relating the 
dune velocity to the length of the dunes. This fit was sur- 
prisingly good and there seems to be a negligible influence 
of the saturation length on dune velocities. The transversal 
and the longitudinal cuts showed that a parabola flts less 
well to numerical and real data than a sin^-function. The 
shape invariance was verified and a characteristic sharp 
edge of the shape at the windward side was observed in 
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Fig. 12. The longitudinal cut through the dune for different 

diffusion constants Cdiff- The length of the slip face and the 
slope at the brink decrease with increasing diffusion constant. 
The shear velocity is u* = 0.5ms~^. 
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Fig. 13. The transversal cut through the windward side of the 
dune for different diffusion constants Cdiff - The height of the 
dune decreases with increasing diffusion constant. The shear 
velocity is w* = 0.5ms~^. 



the calculations. The model was applied to the same dune 
for different diffusion constants. The diffusion constant 
strongly influences tlie form of the final shape. Finally 
it was shown that barchanoids can be obtained from lon- 
gitudinally and laterally connected barchans. 
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Fig. 14. A simulation with quasi-periodic boundary condi- 
tions. There is too much sand to build isolated barchans. The 
barchans are connected in the longitudinal and lateral direc- 
tion. The shear velocity is w* = 0.5ms~'^ 
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Fig. 15. Later state of the same simulation with quasi-periodic 
boundary conditions. Now the barchans are connected only in 
the lateral direction. The shear velocity is w* = 0.5ms~^ 
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nal and lateral direction. 



7. S. Hastenrath. The barchan dunes of southern Peru revis- 
ited. Zeitschrift fiir Geomorphologie, 31-2:167-178, 1987. 

8. H. J. Herrmann and G. Sauermann. The shape of dunes. 
Physica A, 283:24-30, 2000. 

9. P. Hersen, S. Douady, and B. Andrcotti. Relevant length- 
scale of barchan dunes. Phys. Rev. L., 89:264301, 2002. 

10. P. A. Hesp and K. Hastings. Width, height and slope 
relationships and aerodynamic maintenance of barchans. 
Geomorphology, 22:193-204, 1998. 

11. A. D. Howard and J. B. Morton. Sand transport model 
of barchan dune equilibrium. Sedimentology, 25:307-338, 
1978. 

12. D. Jakel. Die Bildung von Barchanen in Faya- 
Largeau/Rep. du Tchad. Zeitschrift fiir Geomorphologie 
N.F., 24:141-159, 1980. 

13. Jose A. Jimenez, Luis P. Maia, Jordi Serra, and Jader 
Morals. Aeolian dune migration along the Ceara coast, 
north-eastern brazil. Sedimentology, 46:689-701, 1999. 

14. G. Kocurek, M. Townsley, E. Yeh, K. Havholm, and M. L. 
Sweet. Dune and dune- field development on Pardre Is- 
land, Texas, with implocations for interdune deposition 
and water-table-controlled accumulation. Journal of Sed- 
imentary Petrology, 62-4:622 635, 1992. 

15. K. Kroy, G. Sauermann, and H. J. Herrmann. Minimal 
model for sand dunes. Phys. Rev. L., 68:54301, 2002. 

16. K. Lettau and H. Lettau. Bulk transport of sand by 
the barchans of the Pampa de La Joya in southern Peru. 
Zeitschrift fiir Geomorphologie N.F., 13-2:182-195, 1969. 

17. K. Lettau and H. Lettau. Experimental and micromete- 
orological field studies of dune migration. In H. H. Let- 
tau and K. Lettau, editors. Exploring the world's driest 
climate, pages 110-147. Madison, Center for Climatic Re- 
search, Univ. Wisconsin, 1978. 

18. H. Momiji and A. Warren. Relations of sand trapping 
efficiency and migration speed of transverse dunes to wind 
velocity. Earth Surface Processes and Landforms, 25:1069- 
1084, 2000. 

19. H. Nishimori, M. Yamasaki, and K. H. Andersen. A simple 
model for the various pattern dynamics of dunes. Int. J. 
of Modem Physics B, 12:257-272, 1999. 



20. M. Sarnthein and E. Walger. Der aolischc Sandstrom aus 
der W-Sahara zur Atlantikkiiste. Geologische Rundschau, 
63:1065-1087, 1974. 

21. G. Sauermann. Modeling of Wind Blown Sand and Desert 
Dunes. PhD thesis, University of Stuttgart, 2001. 

22. G. Sauermann, J. S. Andradc, L.P. Maia, U.M.S. Costa, 
A.D. Araiijo, and H. J. Herrmann. Wind velocity and sand 
transport on a barchan dune. Geomorphology, 2003. 

23. G. Sauermann, K. Kroy, and H. J. Herrmann. A continuum 
saltation model for sand dunes. Phys. Rev. E, 64:31305, 
2001. 

24. G. Sauermann, P. Rognon, A. Poliakov, and H. J. Her- 
rmann. The shape of the barchan dunes of southern Mo- 
rocco. Geomorphology, 36:47-62, 2000. 

25. V. Schwammle and H.J. Herrmann. Modelling transverse 
dunes, cond-mat/0301589, 2003. 

26. M. C. Slattery. Barchan migration on the Kuiseb river 
delta, Namibia. South African Geographical Journal, 72:5- 
10, 1990. 

27. M. S0rensen. An analytic model of wind-blown sand trans- 
port. Acta Mechanica (SuppL), 1:67-81, 1991. 

28. J. M. T. Stam. On the modelling of two-dimensional aco- 
Uan dunes. Sedimentology, 44:127-141, 1997. 

29. J. H. van Boxel, S. M. Arens, and P. M. van Dijk. AeoHan 
processes across transverse dunes i: Modelling the air flow. 
Earth Surf. Process. Landforms, 24:255-270, 1999. 

30. P. M. van Dijk, S. M. Arens, and J. H. van Boxel. Ae- 
olian processes across transverse dunes ii: Modelling the 
sediment transport and profile development. Earth Surf. 
Process. Landforms, 24:319-333, 1999. 

31. I. J. Walker. Secondary airflow and sediment transport in 
the lee of a reversing dune. Earth Surface Processes and 
Landforms, 24:437 448, 1998. 

32. W. S. Weng, J. C. R. Hunt, D. J. Carruthers, A. Warren, 
G. F. S. Wiggs, I. Livingstone, and I. Castro. Air flow and 
sand transport over sand-dunes. Acta Mechanica (SuppL), 
2:1-22, 1991. 

33. G. F. S. Wiggs, I. Livingstone, and A. Warren. The role 
of streamline curvature in sand dune dynamics: evidence 
from field and wind tunnel measurements. Geomorphology, 
17:29-46, 1996. 

34. F. K. Wippermann and G. Gross. The wind-induced shap- 
ing and migration of an isolated dune: A numerical exper- 
iment. Boundary Layer Meteorology, 36:319-334, 1986. 

35. O. Zeman and N. O. Jensen. Progress report on modeling 
permanent form sand dunes. Ris0 National Laboratory, 
M-2738, 1988. 



